Modeling dielectric half-wave plates for CMB polarimetry using a Mueller matrix 

formalism 
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We derive an analytic formula using the Mueller matrix formalism that parameterizes the non- 
idealities of a half-wave plate (HWP) made from dielectric AR-coated birefringent slabs. This 
model accounts for frequency-dependent effects at normal incidence, including effects driven by the 
reflections at dielectric boundaries. The model also may be used to guide the characterization of 
an instrument that uses a HWP. We discuss the coupling of a HWP to different source spectra, 
and the potential impact of that effect on foreground removal for the SPIDER CMB experiment. 
We also describe a way to use this model in a map-making algorithm that fully corrects for HWP 
non-idealities. 
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1. Introduction 

CMB polarization encodes information about Infla- 
tion [l|, re-ionization 0, and the large-scale structure 
of the universe 0] . Manyupcoming experiments includ- 
ing SPIDER 0, EBEXll, POLARBEAR Q, Keck 0, 
ABS [1] and others will use a half-wave plate (HWP) 
to modulate the polarization state of light from the sky 
before measurement by the detectors. Such HWPs may 
be periodically stepped (eg between maps of the same 
part of the sky) to reduce the effect of beam asymme- 
tries and instrumental polarization of optical elements 
skyward of the HWP. Alternatively, they can be contin- 
uously rotated to modulate the signal and to also reject 
atmospheric variations and 1 /f noise . 

An ideal HWP with one of its crystal axes oriented 
at an angle 0hwp to the plane of polarization of inci- 
dent light rotates that polarization plane by 2dhwp as 
the light passes through it. Real HWPs made from bire- 
fringent materials have several important non-idealities. 
Since the phase delay is only a half- wave at a single fre- 
quency, the exact angle through which a HWP rotates 
the polarization state is frequency-dependent. Addition- 
ally, reflections from the material interfaces reduce trans- 
mission and induce non-ideal rotation; to minimize these 
effects anti-reflection (AR) coatings are typically used on 
the surfaces of the HWP. However, AR coatings are fre- 
quency dependent and do not fully eliminate these non- 
ideal behaviors. 

These effects have been treated in a variety of ways by 
other authors. O'Dea et. al. [l^ and Brown et. al. [ll| 
used a Mueller matrix formalism to parameterize a HWP 
and polarized detector, but did not connect their param- 
eterization with a physical model of a HWP. Savini et. 
al. [13 described in detail a physical model for stacks 



of dielectric birefringent materials. Their model works 
directly with the electric fields, which allows it to han- 
dle the input and output polarization state of the light, 
multiple reflections from dielectric interfaces, and the fi- 
nite bandwidth over which the stack is a half-wave re- 
tarder. Matsumura [l^ modeled a multiple-layer HWP 
using a similar approach. We employ the methods of 
Savini et. al. [l^ in this paper, using Jones and Mueller 
matrix methods to derive exact couplings from sources 
(of arbitrary but known spectra) on the sky to a de- 
tector, through a non-ideal HWP made of any number 
of dielectric layers. Our resulting model can be evalu- 
ated quickly in observing simulations such as those done 
for SPIDER [l^ . an important property for upcoming 
CMB experiments. Additionally, we show how to use this 
model to correct for known HWP non-idealities during 
the mapmaking process (using an input data timestream 
to make a map of intensity and polarization on the sky) , 
greatly reducing the potential impact of systematic ef- 
fects induced by such non-idealities. 

2. The HWP Mueller Matrix 

Our goal is to calculate the combined Mueller matrix 
of a HWP plus polarized-detector system for arbitrary 
orientations of the instrument and HWP relative to the 
coordinates defining the Stokes parameters of the incom- 
ing radiation, including the effects of reflections at the 
various dielectric interfaces and the effect of band aver- 
aging. A partially polarized detector can be modeled as 
a partial linear polarizer with Mueller matrix Mpo; fol- 
lowed by a total power detector that is sensitive only 
to the / Stokes parameter. The total Mueller matrix M 
is found by combining the Mueller matrices of the par- 
tial polarizer and the HWP, taking rotational alignments 
into account. The total detected power d for incident 
light with a Stokes vector S = [/, Q, U, V] is given by 
summing over the top row the product of M and S [Tsj , 
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d = I Mil + QMiQ + UMiu + VMiv, 



(1) 
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where the Mxx are elements of the top row of the total 
Mueller matrix. 

Since physical models of the action of a bircfringent 
HWP are constructed using electric fields, we start the 
calculation in the Jones formalism. A Jones matrix is a 
2-by-2 matrix of complex numbers that that describes the 
action of an optical system on the x- and y-components 
of the electric field of an incident plane wave. 

We start by considering the Jones matrix of a general 
retarder, 



iretif) 
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where a(/), &(/), and Mf) are real, and ei(/) and e2(/) 
are small and complex To convert this to a Mueller 
matrix, we follow Jones et. al. [Tsl and use 



Mij — -trace(cri JcTj J^) 



(3) 



from Born and Wolf [l6[, where Ci are the Pauli matrices. 
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This will yield the Mueller matrix as a function of fre- 
quency Mr et{f) of the HWP. Since Mueller matrices 
can be band-averaged, we integrate this single-frequency 
Mueller matrix against a CMB or foreground spectrum 
S{f), as well as the detector passband F{f). This gives 
the band-averaged Mueller matrix 
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A. Single-Plate HWP 

For a HWP made from a single layer of bircfrin- 
gent dielectric material, x- and y-polarized states de- 
fined in the crystal axes cannot couple into each other. 
This means that in the HWP Jones matrix of Eq. [51 
ei(/) = £2(7) = 0. This leaves only three parame- 
ters a{f), b{f) and </>(/) that are necessary to completely 
characterize the HWP. The transmission coefficients a(/) 
and vary with frequency because of the passband of 
the AR coating and the interference of multiple reflec- 
tions inside the birefringent layer. The relative phase 



delay 4){f) also varies with frequency because the path 
length difference for polarization states traveling along 
the slow and fast crystal axes is (n^ — nf)d ~ ^^j- 

Applying Eq. [3] to the Jones matrix of a single-plate 
HWP 



Jret(/) = 

jives the corresponding 



a{.f) 
Mueller matrix Mret(/) 



i(a2 + fe2) i(a2_52) 

|(a2-62) |(a2 + 62) 









abcos{(j)) —absm{(j)) 

absm((f)) abcos((j)) 



(6) 



(7) 



where a, 6, and are all functions of frequency. This 
reduces to the result given in Tinbergen [l7| for an ideal 
retarder (a = 6 = 1). 

In general, the band-averaged HWP Mueller matrix 
Mhwp calculated using Eq. [5] will not be of the same 
form as Mret(/) and will have four (rather than three) 
independent non-zero elements. 
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For the specific case of the SPIDER 145 GHz HWP, we 
have calculated this Mueller matrix as a function of fre- 
quency, and plotted it in Fig. [TJ Since the CMB is not 
expected to be circularly-polarized, in the case where the 
subsequent detector and optical system does not induce 
sensitivity to circular polarization the s parameter will 
not be relevant for CMB polarimetry. 

B. Multiple-Plate HWP 

For a multiple-layer HWP, we cannot simply multiply to- 
gether the Mueller matrices of several single-layer HWPs, 
because this would not account for multiple reflections 
among the interfaces between each layer. We instead re- 
turn to the original Jones matrix of a general retarder in 
Eq. [21 For a multiple-layer HWP, ei and £2 will be small, 
but still non-zero. This means that all seven of the pa- 
rameters of the Jones matrix, a{f), b{f), (j){f), Re[ei(/)], 
Im[ei(/)], Re[e2(/)], and Im[e2(/)] are required to char- 
acterize the HWP. After applying Eq. [31 we obtain the 
corresponding Mueller matrix 
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Fig. 1. Calculated Mueller matrix elements (see Eq. [8]for where these elements lie in the matrix) as a function of 
frequency of an optimized SPIDER 145 GHz HWP (solid blue lines). The dot-dash red lines show the Mueller matrix 
elements of a HWP retarder with a perfect AR coating for comparison. The grey box shows the SPIDER 145 GHz 
band. 
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a' + - |ei|" - |e2| 
b ■ (Rc[ei] cos(^) + Im[ei] sin((/))) 
b ■ (Rc[ei] sin((/)) — Im[ei] cos((^)) 
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This reduces to the single plate result of Eq. [7] if ei = 
£2 — 0. Note that none of the elements of this matrix 
have the same functional form. This means that after 
applying Eq. [5] to band-average against the source and 
detector spectra, all 16 elements of the resulting Mueller 
matrix will be independent. Since the CMB is not ex- 
pected to be circularly-polarized and the subsequent op- 
tical and detector system should have no sensitivity to 
circular polarization, only the top-left 9 elements of the 
Mueller matrix will be relevant for CMB polarimetry. 



3. Rotating the Instrument and HWP 

Jones ct. al. [l^ modeled a polarization-sensitive de- 
tector as a rotatablc instrument with a partial-polarizer 
followed by a total power detector. The Jones matrix of 
a vertical partial polarizer is 



■ipol 
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(10) 



which can be turned into a corresponding Mueller matrix 
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the detected radiation is given by the coupling to / in 
the product Mpoi M^,; for an ideal polarizer with rj = 1 
and 6 = and an arbitrary instrument angle ipmst, the 
detector signal is 

d=^{I + Q cos{2ipi„st) + U am{2ijinst)) . (13) 

The addition of a HWP to the instrument can be mod- 
eled with the Mueller matrix product 



M = M 



poi M_ 



'\hwp Mg M 



(14) 



where is the rotation matrix by the HWP angle fl/i^p, 
and is the rotation matrix by the detector orienta- 
tion angle ^det ■ Fig. [5] illustrates the definition of these 
angles; while only two angles are needed to define the 
single-detector problem, we use three here to make the 
problem more straightforward to visualize, and to easily 
accommodate calculations of focal planes with detectors 
at multiple orientation angles. Here we consider only a 
single detector with ^det = 0. 

The top row of the resulting Mueller matrix for a 
single-plate HWP is 

Mil = ^ [T{t]^ + 6^) +pcos{20h^p)i^^ - 6^)] 

MjQ = T sin{2il;inst) + G cos{2ip,„st) 
Miu = -J'cos(2V'mst) + 0sin(2V'mst) 



Miv = -ssin(26l,™p)(?7^ -,5^), 
where T and Q are 

- ^psm{2dhu,p){if +3"^) 
Q = i [T + c + (T - c) cos(40,,„p)] (772 - <52) 



(15) 



(16) 



(17) 



This reduces to the result given in Jones et. al. [l5| in 
the limit of no HWP (T = c= l, p=s = 0), and in 
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Fig. 2. Our definitions of the detector, HWP, and instru- 
ment angles. The instrument angle is defined relative to 
a fixed reference on the sky that determines the abso- 
lute orientations of Q and /7, and the detector and HWP 
angles are defined relative to the instrument. 



the case of an ideal HWP (T = 1, c = -1, /? = s = 0) 
reduces to 
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(18) 



Armed with the band-averaged Mueller matrix repre- 
sentations of the HWP plus detector system given by 
Eq.El we can use Eq. [T]to calculate the detector output 
as a function of input Stokes parameters. 

An example of the output from our model, along with 
the percent-level differences between that and a naiive 
treatment, is shown in Fig. [3] for the specific case of 
the SPIDER experiment, which uses a single birefringent 
sapphire HWP with quarter-wave AR coats. 

4. Detector Pair Summing and Differencing 

Many CMB polarization experiments use pairs of 
Polarization-Sensitive Bolometers (PSBs) located at the 
same point on the focal plane [l^. Both detectors in a 
pair view the same patch of the sky, but detector B is 
oriented at 90° with respect to detector A. Detector dif- 
ferencing within a pair reduces common-mode noise while 
retaining sensitivity to linear polarization. The sum of 
a pair measures the / Stokes parameter. We take Eq. [1] 
with M given by Eq.[T3]with ^^et = as a model for the A 
detector timcstream df. The B detector timestream df 
can be similarly calculated by setting S^det — 90° ; we note 
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Detector Signal: Mueller Matrix Method (ideal HWP) - (Simulated Real HWP) 




Instrument Angle [deg] Instrument Angle [deg] 

Fig. 3. Simulated detector output for the SPIDER instrument as a function of HWP and instrument angle for a Q- 
polarized CMB dB/dT source. The left panel shows the detector output calculated using our band-averaged Mueller 
matrix formalism. The right panel shows the percent-level differences between an ideal HWP and the simulated 
non-ideal HWP. 



however that the Jones matrix of a horizontal polarizer, where F' and Q' are 



1^ = 

-'pol 
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(19) 



is related to that of a vertical polarizer (Eq. [T0| via the 
substitutions rj ^ 5 and S ^ rj. We can thus just make 
those substitutions in Equations fTSl through [iTl to find a 
model for the B detector timestream. We then construct 
the sum and difference timestreams 



dt+df, 
df-df. 



(20) 



Only ((5^ + rf') terms will remain in the sum timestream 
and only ((5^ — if') terms will remain in the difference 
timestream. The matrix elements for the sum timestream 
with a single-plate HWP are therefore 



0. (21) 



Even with the non-idealities of the HWP, the sum 
timestream coupling to intensity is independent of HWP 
angle. There is also a small coupling to linear polariza- 
tion. The matrix elements for the difference timestream 
are 

Mfg" = [T' sin(2V'm,,t) + g' cos(2^,„,t)] iv^ ~ 
Mfjf = cos(2V'.„..t) + g' sin(2^„,t)] - 5^) 



M 
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IV 



s sm{2d hu,p)i'n^ - S^), 



(22) 



T' = --(T-c)sin(40,,^p) 



g' 



-[T + c+iT-c) cos(40,,^p)] . 



(23) 
(24) 



The difference timestream unfortunately has small cou- 
plings to intensity and circular polarization. Note that 
for both the sum and difference timestreams, detector 
cross-polarization only shows up as an overall factor of 
(77^ ± 5^), and will not result in leakage between the 
estimates of Q and U. 

5. Application to a Specific Instrument 

The model we derived requires input values for the HWP 
parameters a(/), 6(/), </>(/), ei(/), and e2(/). We can 
calculate these parameters from first principles, or mea- 
sure them using a polarized Fourier transform spectrom- 
eter (FTS) with the assembled instrument. 
A. Physical Modeling 

To calculate the HWP parameters from first principles, 
we use a physical optics model similar to the one in Savini 
et al. The model extends the 2-by-2 matrix for- 

malism described by Hecht and Zajac [1^ for modeling 
multiple dielectric layers of isotropic materials to a 4- 
by-4 matrix formalism for multiple layers of potentially 
birefringent materials. The model uses the electromag- 
netic boundary conditions at the interface between each 
layer of material to map the incident electric and mag- 
netic fields onto the transmitted fields. This fully treats 
multiple reflections and interference effects, and can also 
handle lossy materials. 

We use this model to calculate the elements of the 
HWP Jones matrix shown in Eq. [5] by calculating the 
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transmitted electric field amplitude Eout at a frequency / 
with an incident electric field amplitude Ein, for both x- 
and y-polarized incident waves. We assemble the trans- 
mitted amplitudes into the Jones matrix for that HWP, 



«(/) 



Eout 
- Eout 



(25) 



B. Instrument Calibration 

The parameterization we developed can also be used to 
characterize the HWP, with the aim of better under- 
standing the coupling to intensity and linear polariza- 
tion signals from the CMB and foregrounds with other 
spectra. This characterization process depends on the 
manner in which the HWP will be used, ie whether it 
will be continuously rotated or occasionally stepped dur- 
ing the observation process. For both the stepped and 
continuously rotating cases, a polarized FTS can be used 



to measure the frequency and polarization response of 
the instrument and HWP. Taking the Fourier transform 
of the detector timestream as the FTS mirror moves at 
a constant speed will yield the spectral response of the 
combined instrument |2C|. 



1. Stepped HWP 

As an illustration of how the calibration process of an 
instrument with a stepped single-plate HWP might work, 
consider the case where we align the instrument with the 
FTS such that ^det = V'msi = 0°. The FTS spectrum 
S{Ohwp,f ) at a particular HWP angle can be modeled 
using Eq. [1] with the Mueller matrix elements given in 
Eq.[T4l Before proceeding, it is useful to change variables 
from T] and 6 to the optical efficiency and the polarization 
efficiency 7. Here we define the optical efficiency as 



i(77^ + (S^) and follow Jones et. al. 



polarization efficiency 7 



15| and take the 
Setting ijj.nst = and 



collecting terms in Eq. [T] gives 



5(^W,/) = (^/+^gcos(40;,„,p)^T(/)e(/)+ QQ + i[/sin(40^„p)^r(/)e(/)7(/) 
iQcos(20^.„p) + ^Usm{2ehu,p)] p{f)e{f) + I cos{2ehu.p)p{f)e{fh{f) 



+ i [Q(l - cos(40^^p)) - C/sin(40„^p)] c{f)e{fh{f), 



(26) 



r 



where /, Q, and U are the Stokes parameters of the light 
coming out of the FTS. Placing a polarizing grid at the 
output of the FTS produces light with 1=1 and Q = 1. 
We call spectra taken in this configuration Sq. Rotating 
the grid by 45° produces light with / = 1 and U = 1; 
spectra taken with this configuration we call 1S45 . 

We have five parameters to constrain (for a single-plate 
system; a multi-plate HWP has more) , so we need at least 



five FTS spectra with different instrument and HWP an- 
gle combinations. Considering noiseless data, we need 
spectra with the polarizing grid at both and 45° to in- 
dependently estimate all of the parameter combinations. 
With this data, we can set up a system of linear equations 
to estimate the HWP and detector parameters. 



( 



?Af+2, /) 



2 + cos(46'i) 1 cos(26ii) 
2 + cos(46'2) 1 cos(26'2) 
2 + cos(46i3) 1 cos(203) 



2cos(26li) l-cos(46'i) 
2cos(26'2) l-cos(46i2) 
2cos(26'3) l-cos(46i3) 



2 sin(46'jv+i) sin(26'jv+i) 2 cos(26'7v+i) - sin(46'Ar_|_i) 

2 sin(46'7v+2) sin(26'7v+2) 2 cos(26'7v+2) - sin(46'Ar+2) 



nmf) 

p{f)e{f) 

pifHfhif) 

c(/)£(/)7(/) 



(27) 



Operationally, estimating this combination of the pa- (detectors -I- HWP) for polarized CMB mapmaking. 
rameters is sufficient to characterize the full instrument 
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However, since the detector optical efficiency e{f) ap- 
pears as an overall factor in these equations, we cannot 
independently estimate the HWP-only parameters T{f), 
p{f), and c(/). To check the HWP physical optics model, 
we would need separate FTS spectra without the HWP 
in place to estimate the detector optical efficiency. 

In this illustration, we do not consider possible pa- 
rameter covariances or biases induced by detector noise 
and non-optimal instrument and HWP angles. One pos- 
sible way to handle this would be to use FTS spectra 
taken at many HWP and instrument angles to obtain 
least-squares estimates of the calibration parameters and 
errors as a function of frequency. This would also allow 
a test in the sense on whether the performance of the 
HWP and instrument is well-described by our formalism. 
2. Rapidly-Rotating HWP 

If an instrument designed for a rapidly-rotating single- 
plate HWP has a mode in which the HWP can be 
stepped, then obviously the characterization procedure 
in the previous section may be used to estimate the HWP 
parameters. 

As a substitute or complement, the HWP can be char- 
acterized in the instrument as it is rapidly rotating. For 
a single-plate HWP, taking ipinst = 0° and collecting the 
constant, 2fhwp and 4fhwp terms in Eq.[T] gives 



d{6hwp) 



Tel + -{T + c)e-fQ 



-{T-c)ejU 



pe-fl + -peQ 



sin46',,„p - 
cos 26hwp 



-{T^c)ejQ 



se'yV + ^P^^ 



(28) 

cos 4,9hnjp 

sin 26hwp- 



This shows that only polarization information is con- 
tained in the '^Jhwp component of the detector 
timestrcam. A lock-in at -ifawp on the timestream di 
will produce companion Xi and timcstreams, where 
Xi is the sine wave lock-in output timestrcam, and yi is 
the cosine timestrcam. The Q and U timcstreams in in- 
strument coordinates (since ^det = V'mst = 0°) may be 
estimated using 



Q^ = 

U, = 



(T - c)£7 
2 



(29) 



This means that only one calibration factor is necessary 
to make polarization maps with a rapidly-rotating single- 
plate HWP. Doing this as a function of frequency with an 
FTS allows for the estimation of this calibration factor 
for different source spectral shapes. 

6. Polarized Mapmaking 

In many CMB polarization experiments, the instrument 
continuously scans the sky, and the detector timestreams 
are used to estimate maps of the /, Q, and U Stokes pa- 
rameters. Since the CMB is not expected to be circularly- 
polarized, we take V = Q. To estimate the maps, a 



mapmaking algorithm such as the one in Jones et. al. 
jl5| assumes that without cross-polarization, the detec- 
tor timestream d, can be modeled as 



= i (/ + Qcos(2V'i) 



U sin(27/;i)) + 



(30) 



where ipi is the instrument angle timestrcam, and rii is 
atmospheric and detector noise. The algorithm then 
combines this detector timestream with the pointing 
timestream of the telescope, and uses an iterative ma- 
trix method to compute maximum likelihood maps of /, 
Q, and U , as well as an estimate of the noise timestream. 
The algorithm described in Jones et. al. [l5| also can 
treat cross polarization in the detector. 

As an illustration of the principle of the method ap- 
plied to mapmaking with a HWP, we consider TV noiseless 
detector samples taken at different instrument angles ipi 
and HWP angles Ot when the telescope was pointed at 
the same spot on the sky. To estimate the Stokes pa- 
rameters at that point on the sky, we set up a system of 
linear equations using the detector timestream model. 



/ 

Q 

u 



(31) 



Here c? is a A'^-clement column vector of detector sam- 
ples, and A is a A^-by-3 matrix containing the detector 
and HWP model. If we take Eq. [30] and model an ideal 
polarized detector with an ideal HWP, the A matrix is 



1 cos(2(V'i -^26*1)) sin(2(^/'i +26I1)) 
1 cos(2(V^2 + 26I2)) sin(2(V'2 + 26I2)) 

1 cos(2(V'w -l-26ijv)) sin(2(i/'Ar + 26'A,)) 



.(32) 

To estimate the polarization state of the incident light, 
we invert the A matrix. 



/ 

Q 

u 



A-^d. 



(33) 



To treat the effects of the HWP non-idealities, we can 
use the timestrcam model of Eq.[T]and the Mueller matrix 
elements from Eq. [T3] to set up the relevant system of 
linear equations. In this case the A matrix becomes 



M„(0i,^i) 



M/q(0i,Vi) 

M/q(02,V'2) 



Miu {02,^^2) 



Mii{dN,i^N) MlQi9N,i'N) Miu{0n,i1^n) 

(34) 

Even though this matrix has more complicated elements, 
Eq. [33] still holds and we can still estimate the polar- 
ization state of the incident light with a matrix inver- 
sion. Thus, simply changing the functional form of the A 
matrix in the mapmaking algorithm allows it to exactly 
compensate for the HWP non-idealities in our model. 
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7. Application to the SPIDER HWP 

A. Mueller Matrix 

The HWPs in the SPIDER mstrument will consist of a 
single 330 mm diameter single-crystal sapphire plate for 
each telescope with a quarter-wave quartz layer attached 
to each sides. Each HWP will be cooled to ^ 4 K. Here 
we consider the prototype HWPs for the SPIDER 145 
GHz band. 

We have measured the indices of refraction of sapphire 
at 5 K to be 3.336 ± .003 and n/ = 3.019 ± .003, 

from 100 GHz to 240 GHz [2l|. Our fused quartz AR 
coats have a room temperature index of refraction of 
riar = 1.951 at 245 GHz [H. We use Eq. [25]to calcu- 
late the Mueller matrix of a thickness-optimized SPIDER 
145 GHz HWP as a function of frequency; the results are 
shown in Fig. [TJ We assumed all materials were lossless, 
though that assumption can easily be relaxed and the 
output used in the Mueller formalism described above 
without modification. 

We then use Eq. [5]to calculate the band-averaged ma- 
trix elements of the SPIDER HWP, assuming a top-hat 
detector spectrum from 130 GHz to 160 GHz. For the 
source spectra we use 

' 1 Flat 
^(/,2.725 K) CMB [23] 

^..^ l./'-''S(/,9.6 K) 

I +0.0935/2-^B(/,16.2 K) Dust [24] 
/^^ Synchrotron [25] 

J- 14 Free-free [26], 

(35) 

where B{f,T) is the blackbody function, as estimates of 
the CMB and astrophysical foregrounds. The results are 
shown in Table [T] The c parameter deviates from ideal- 
ity by almost 5%, which is a relatively large effect. This 
pushes us towards handling the non-idealities through 
calibration and a modified mapmaker as discussed in Sec- 
tion [51 

In addition to being large amplitude, the s parameter 
depends significantly on which source spectrum is used; 
this is not surprising given the form of s across the band, 
shown in Fig. [1] which shows that tailoring s to be near 
zero requires near symmetric placement and weighting 
across the band, since s is large but asymmetric. The 
variations of T, p, and c between the CMB and fore- 
ground sources are not as worrisome. The parameters 
vary among the different sources at less than the 10^"* 
level. 

For comparison, SPIDER is targeting B-modes at the 
r ~ .03 level, and we will remove dust foregrounds by 
combining maps from several frequencies. To reach our 
target i?-mode sensitivity, our goal for dust foreground 
subtraction is to reduce it by ^ 95%. The HWP response 
varies among the different sources at the 10"^ level, 
which suggests that residuals from foreground subtrac- 
tion should be dominated by incomplete dust removal. 





T 


P 


c 


s 


Flat 


0.97389 


0.01069 


-0.95578 


0.00170 


CMB 


0.97396 


0.01069 


-0.95591 


-0.00952 


Dust 


0.97391 


0.01080 


-0.95563 


-0.03598 


Synchrotron 


0.97382 


0.01070 


-0.95565 


0.01280 


Free- free 


0.97388 


0.01069 


-0.95577 


0.00325 


(Ideal HWP) 


1 





-1 






Table 1. Calculated Mueller matrix elements for an op- 
timized cryogenic sapphire HWP with a quartz AR-coat 
for the Spider 145 GHz HWP. We optimized the HWP 
thickness based on our measured cold indices, and the 
AR-coat thickness based on the index in Lamb (2^ . 
The first row shows the HWP parameters averaged 
within the Spider 145 GHz passband. The CMB, Dust, 
Synchrotron, and Free-free rows all are band-averaged 
against the source spectra within the passband. The last 
row shows the parameter values of an ideal HWP for 
comparison. 



not source-varying HWP systematics. Next-generation 
experiments targeting r < .01 will require better fore- 
ground removal, and may be impacted by this effect. 

B. Mapmaking 

A full time-domain simulation of a given instrument's 
scan strategy is necessary to precisely quantify the ef- 
fects of using a conventional mapmaking algorithm with 
a non-ideal HWP, and to see if the modified mapmaker 
described above improves the results. To get a sense for 
what the results might be, we performed a simpler sim- 
ulation. We input a Stokes vector of [100 ^K, 1 /iK, 0, 0] 
to the detector model in Eq. [TJ This Stokes vector is 
representative of the amplitudes of temperature and E- 
mode CMB fluctuations at the angular scales relevant 
for SPIDER. We fixed ^det = i'mst = 0° and gener- 
ated one detector sample at each of the HWP angles 
0°, 1 X Yg3f,2 X jg^Y , . . . , ^mar to make a total of 10 
samples evenly spaced in angle from 0° to 9max- We 
used these simulated detector samples to estimate the Q 
and U Stokes parameters using Eq. [331 We simulated 
the effect of using a conventional mapmaker by using the 
matrix AideaU ^.nd simulated a modified mapmaker by 
using Areai- Wc repeated this process with varying 9max 
to generate Fig. 31 which is a plot of reconstructed Q and 
[/ as a function of Omax- We generated the second pair 
of plots shown in Fig. [4l with an input Stokes vector of 
[100 AiK,0,l AiK,0]. 

Not shown in Fig. [4l is the fact that the conventional 
mapmaker also mis-estimates the / Stokes parameter 
when using samples from a detector with a non-ideal 
HWP. However, most CMB instruments (including SPI- 
DER) will calibrate by cross-correlating the intensity 
maps with the WMAP or upcoming Planck satel- 
lite results, which would change the overall calibration 
to compensate for the / mis-estimation. To approximate 
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lnputl = 100nK, Q=1 nK, U = 0|iK Input I = 1 00 nK, Q = piK, U = 1 nK 




90 180 270 360 '0 90 180 270 360 



0.8 
0.6 

D 

0.2 


-0.2 



90 180 270 

Maximum HWP Angle [deg] 



360 



2 
1.8 
1-6 

D 

1.2 
1 

0.8 





• Standard Mapmaker 
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A 


•*•» 
• • 



90 180 270 

Maximum HWP Angle [deg] 



360 



Fig. 4. Q and I] Stokes parameter estimates made with a conventional mapmaker and a mapmaker modified to handle 
the SPIDER HWP non-idealities. Each estimate was made with simulated detector samples from 10 HWP angles 
evenly spaced from 0° to Omax, where 9max is the horizontal axis of each plot. The conventional mapmaker, plotted 
in dots, misestimates the Q and U Stokes parameters by up to 1 yuK. The main reason for this is the temperature- 
polarization leakage in the HWP. The modified mapmaker, plotted as crosses, completely handles the non-idealities 
modeled in our formalism. 



the effect of this calibration process, we scaled our esti- 
mates of /, Q and U made with the conventional map- 
maker by a factor to force the / estimate to be the same 
as the input value. This was not necessary for the mod- 
ified mapmaker since it accurately estimates the Stokes 
/ term. The reason the conventional mapmaker gener- 
ates the large residuals shown in Fig. |4] is that it does 
not handle the temperature-polarization leakage gener- 
ated by the HWP. 

Quantifying the systematics from using a non-ideal 
HWP with a conventional mapmaker will require full 
time-domain instrument simulations. These will be pre- 
sented for the case of SPIDER in an upcoming paper 




8. Conclusions 

Wc have derived an analytic model for birefringent, AR- 
coated HWPs that simulates the response of an instru- 
ment to different source spectra on the sky at arbitrary 
instrument and HWP angles. We also modeled the sum 
and difference timestreams of a PSB pair looking through 
a HWP, and showed that it should be possible to char- 
acterize the HWP and detectors with lab testing of a 
completed instrument. Motivated by the presence of po- 



tentially significant non-idealities in the HWP, we pre- 
sented a mapmaking algorithm that accounts for these 
known non-idealities of the HWP. This model will be in- 
tegrated into an end-to-end simulation of the SPIDER 
instrument to translate these non-idealities into detailed 
limits on their possible contamination of power spectral 
and parameter estimates. 

After the initial calculation of band-averaged param- 
eters, our Mueller matrix method is analytic and does 
not require repetitive matrix multiplications or repet- 
itive integration over frequency. Therefore, modeling 
detector timestreams is far faster than repeated use of 
a Jones-formalism code. As an example, in 10 sec- 
onds of computer time in Matlab on a laptop, the di- 
rect Jones method can simulate 300 SPIDER detector 
samples, while the Mueller matrix method can simulate 
24,000, a speedup by a factor of ~ 75. A version of the 
Mueller matrix code written in C for use in the instru- 
ment simulation code for SPIDER can simulate 10 million 
detector samples in 10 seconds on the same machine. 
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